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A solid-state electronic nanodevice is an intrinsically open quantum system, exchanging both en- 
ergy with the host material and carriers with connected reservoirs. Its out-of-equilibrium behavior 
is determined by a non-trivial interplay between electronic dissipation/decoherence induced by in- 
elastic processes within the device, and the coupling of the latter to metallic electrodes. We propose 
a unified microscopic description that accounts for both these aspects, enabling to predict various 
steady-state as well as ultrafast nonequilibrium phenomena, nowadays experimentally accessible. 
More specifically, we derive a generalized density-matrix equation, particularly suitable for the de- 
sign and optimization of a wide class of electronic and optoelectronic quantum devices. The power 
and flexibility of this approach is demonstrated with the application to a photoexcited triple-barrier 
nanodevice. 
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The never-ending technological progress in solid-state 
meso- and nanosystems [Tj , such as semiconductor super- 
lattices [2] and quantum-dot structures [3j, has enabled 
to observe several phase-coherence phenomena [U [5J, 
where the wavelike properties of electrons become ap- 
parent. In any realistic quantum device [5J [7J, how- 
ever, phase coherence is hindered by the inelastic scat- 
tering that electrons experience within the host ma- 
terial (phonons, photons, plasmons, etc.) as well as 
by the coupling to external reservoirs; such phenom- 
ena may cause energy dissipation, decoherence, and car- 
rier thermalization, thereby affecting the nanodevice be- 
havior. Within such general framework, two different 
and well-established fields may be identified. On the 
one hand, the progressive reduction of the space- and 
time-scales of new-generation electronic and optoelec- 
tronic devices [7J forces one to replace the traditional 
Boltzmann picture [5] with genuine quantum approaches 
based, e.g., on the density-matrix formalism [2 [7], on the 
Green's function theory [9j, and on the Wigner- function 
picture [10]. On the other hand, a realistic descrip- 
tion of new-generation quantum-transport devices [B] re- 
quires to extend the Landauer-Biittiker treatment of the 
mesoscopic-transport regime [llj . in order to include 
energy-dissipation and decoherence phenomena. Aim of 
this Letter is to provide a unified microscopic framework, 
based on the density-matrix formalism, able to describe 
both the above-mentioned fields. 

The crucial interplay between phase coherence and 
dissipation/decoherence phenomena is often described 
through the equation of motion for the electronic single- 
particle density matrix [5] [7] 
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Here, the first term on the r.h.s. accounts for the 
coherent evolution, dictated by the non-interacting 
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FIG. 1. (Color online) A nanodevice as an intrinsically open 
quantum system, which can exchange energy with the host 
material (environment), and carriers with external reservoirs. 
Its nonequilibrium properties are determined by the non triv- 
ial interplay between these two types of interactions. 



single-particle Hamiltonian H° = "Yl, a \o()e a {a\, where 
| a) denote the single-particle eigenstates (corresponding 
to the energy levels e a ), typically given by the scattering 
states related to the device potential profile. The second 
term encodes dissipation and decoherence processes, 
arising from the energy exchange between the carriers 
and the host material; henceforth we shall refer to 
such term as the carrier-environment (env) coupling. 
Equation (TlJ) applies to a broad variety of problems, 
a remarkable example being the semiconductor Bloch 
equations [7j. 

Importantly, the degree of accuracy of Eq. (IT]) is closely 
related to an appropriate choice of its last term. Indeed, 
oversimplified phenomenological treatments can lead, for 
instance, to a violation of the positive-definite character 
of the density-matrix operator p = J2 a a \ a i)Pa 1 a 2 ( a 2\, 
which is a mandatory prerequisite of any quantum- 
mechanical time evolution. To this end, a general pre- 
scription is to express various energy-relaxation coupling 
mechanisms via suitable Lindblad superoperators [12] ; In 
this way, the last term in (flj can be written in operatorial 



form as 



dp 
dt 
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where A s denotes the Lindblad superoperator related 
to the s-th interaction mechanism with the host material. 

Besides the environment, an electronic/optoelectronic 
nanodevice is also coupled to so-called metallic elec- 
trodes, i.e., to macroscopic charge reservoirs. This makes 
any realistic nanodevice an intrinsically open quantum 
system, exchanging energy with the environment and car- 
riers with various reservoirs, as schematically depicted 
in Fig. [11 However, these two interaction channels can- 
not be treated on equal footing. Indeed, the coupling to 
reservoirs cannot be accounted for via a term like (pL 
which necessarily implies that d(trp)/dt = 0. Different 
approaches are thus needed to account for the fact that, 
while p remains a positively-defined operator, its trace is 
in general not preserved. In the coherent regime, where 
inelastic coupling with the environment is negligible, the 
Landauer-Buttiker formalism [11] provides a successful 
description of the system steady state in the presence of 
external reservoirs. However, when inelastic scattering 
mechanisms are present, this treatment becomes inade- 
quate. An alternative strategy is the Wigner-function 
approach jlOj . where dissipation/decoherence phenom- 
ena due to the environment are usually described within 
the relaxation-time approximation, whereas the presence 
of the reservoirs is encoded via suitable spatial bound- 
ary conditions. However, such classical-like approach has 
been recently shown |13j to lead to negative carrier prob- 
ability densities. 

A general and physically consistent description, in- 
cluding dissipation and decoherence phenomena due to 
the open character of the electronic system, is therefore 
still lacking. In the present Letter we propose a unified 
and physically reliable framework for the microscopic de- 
scription of open quantum systems, which enables us to 
account for both dissipation/decoherence due to the en- 
ergy exchange with the environment and carrier transfer 
from/into the reservoirs. We shall show that this ap- 
proach is particularly suitable for the design and opti- 
mization of new-generation quantum devices with open 
spatial boundaries, corresponding, e.g., to the case of a 
semiconductor nanodevice inserted into an electric cir- 
cuit, schematically depicted as the N reservoirs in Fig.[T] 

The explicit form of the desired system-reservoir cou- 
pling superoperator should fulfill three basic require- 
ments. It should (i) preserve the positivity of the density- 
matrix operator />; (ii) not induce phase coherence on 
the carrier system, due to the thermal (i.e., fully inco- 
herent) nature of the external reservoirs; (iii) reduce to 
the standard injection-loss structure in the semiclassical 



limit (p ai a 2 = f ai S ai a 2 ) 0, i-e. 

=5„-ir/a = -rr(/a-/r), (3) 
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where T r ° s denotes the inverse timescale for f a to reach 
the steady-state distribution /^ cs in the absence of 
carrier-environment coupling. 

Adopting a fully operatorial notation, the proposed ex- 
tension of the density-matrix equation (nj), compatible 
with all the above requirements, reads 
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is the density-matrix operator encoding the distribution 
of the incoming reservoir electrons, with fy denoting the 
equilibrium or quasicquilibrium carrier distribution of the 
j-th reservoir [14], whereas 



B k3 =\a k] )jT^W\ 



(7) 



are Lindblad-likc superoperators describing the coupling 
between the generic free-particle state \k 3 ) incoming from 
the j-th reservoir and the corresponding single-particle 
scattering states |a#) of the nanodevice. Differently 
from the genuine Lindblad form in ([2]), describing the 
energy exchange with the environment, the first term on 
the r.h.s. of Eq. ^ involves the density-matrix operator 
p° of the external reservoirs. This feature, in sharp 
contrast with the superoperator in S, makes Eq. Q 
inhomogeneous, implying that tip is not conserved, as 
expected in a system where the number of particles may 
change with time. Nevertheless, the positive-definite 
character of p is ensured by the Lindblad-like form of 
the proposed coupling term. It is also worth noticing 
that in our approach the local (i.e., classical- like) 
boundary-condition treatment of the system-reservoir 
interaction employed in the Wigner formulation [13] 
has been replaced by a non-local quantum-mechanical 
coupling. 

It is useful to rewrite the operatorial equation Q 
within the single-particle basis a. By introducing the 
density-matrix operator p rcs = V ■ E feJ K4/a„, («# I ^ 



p°, and by inserting Eqs. ^ and ^ into Eqs.Q-^ one 
obtains 



dp ai 
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where the compact notation a fcJ — > a has been employed. 

In order to illustrate the implications of the system- 
reservoir coupling [second line in Eq. Q], let us start 
by analyzing the mesoscopic regime, i.e. the case where 
inelastic scattering mechanism are negligible and thus 
the environment coupling term in Eq. (ISj) can be omit- 
ted. Then, the steady-state (ss) solution of Eq. ([8| is 
easily obtained as p^ ia2 = p^ a2 = f°5 aiCX2 , and the 
Landauer-Buttiker result [TT] is recovered: the steady- 
state single-particle density matrix p s ^ ia2 of the device is 
diagonal, thus preventing any phase coherence transfer 
between reservoirs and device. Its diagonal values coin- 
cide with the distribution f° of the carriers injected from 
the reservoirs. Notice that such steady-state solution 
Pa 1 a 2 i s independent of the value of the device-reservoir 
coupling constants T™ s . However, our approach goes be- 
yond the steady-state regime and enables us to address 
ultrafast-dynamics phenomena, nowadays accessible via 
time- resolved optical experiments. In particular, if one 
initially "prepares" a device electron in a coherent super- 
position of the single-particle states \a), characterized by 
a device density matrix ~p aia2 at t = 0, its subsequent 
time evolution is given by: 



Paia 2 {t) — P ai a 2 ' \Pa 1 a 2 



Pa\Oi 2 



(9) 



where S Q = e a + ihT T ° s /2 can be regarded to as 
the device single-particle self-energy e Q "dressed" by 
an imaginary life-time contribution induced by the 
device-reservoir interaction. Thus Eq. ^ enables one to 
interpret the coupling T a = T a . as the inverse time- 
scale over which the diagonal density-matrix element 
p aa reaches the steady state value p s ^ a , if only injection 
from the j-th reservoir were present. Furthermore, the 
time evolution in |9| shows the interplay between phase 
coherence (encoded in the non-diagonal density-matrix 
elements) and reservoir-induced dissipation/decoherence 
processes [T5] . 

Let us now include the presence of inelastic interaction 
mechanisms with the environment (phonons, photons, 
plasmons, etc.). In this case, the time evolution is 
given by the generalized density-matrix Eq. pi). As 
shown recently [16 , by adopting an alternative Markov 
procedure, it is possible to perform a fully microscopic 
derivation of the Lindblad superoperators in (pi); the 
latter, written within the basis a, involve non-diagonal 
scattering rates expressed via a generalized Fermi's 



golden rule. The highly non-trivial structure of such 
scattering superoperators does not allow for an analyt- 
ical treatment of the problem; numerical solutions are 
needed and will be presented here below for a realistic 
nanodevice. Nevertheless, some qualitative remarks can 
be made already at this stage. Indeed, the existence of 
such microscopic derivation of the Lindblad superop- 
erators implies that a parameter-free characterization 
of the system-environment is possible for any given 
device. Then, when a nanodevice is disconnected from 
the external circuit, it is possible to determine the 
steady-state density matrix p c ™ a2 reached in the absence 
of external reservoirs, and to qualitatively describe such 
thermalization process via partially phenomenological 
coupling constants r° nv = £ s r™ v > s , where T™ v ' s is 
operatively defined as the inverse timescale that the 
system would take to reach the steady-state solution, if 
only the s-th interaction mechanism with the environ- 
ment were present; the specific value of r° nv,s depends 
strongly on the nanodevice geometry as well as on the 
specific interaction mechanism considered. 

When the system is further connected to the exter- 
nal circuit, the steady-state solution is the result of a 
non-trivial interplay between dissipation versus decoher- 
ence processes induced by energy exchange with the en- 
vironment, and carrier exchange versus thermalization 
induced by the external reservoirs. It is then quite nat- 
ural to identify two limiting regimes: for T™ v 3> r^ cs 
the effect of nonelastic interaction processes is domi- 
nant over the reservoirs one, and in this limit we get 
Pa ia2 = /°™a 2 ' i- e -' tne en ~ ect °f the scattering processes 
within the device is such to maintain the electron gas in 
thermal equilibrium with the host material; in contrast, 
for r™ v <C T r ° s the effect of inelastic scattering is negligi- 
ble, and in this limit the mesoscopic steady-state solution 
Pa 1 a 2 = Paia 2 ^ s recovered, i.e., the latter is essentially 
determined by the external reservoirs. 

In order to show the power and flexibility of the 
proposed quantum-transport formalism, we have ap- 
plied this simulation strategy to the investigation of 
dissipation-versus-decoherence phenomena in a realis- 
tic nanodevice [T7]; to this end, we have considered 
a symmetric GaAs/AlGaAs triple-barrier structure (see 
Fig. [2]), connected to two electrodes, tailored such to al- 
low for a significant interwell tunneling corresponding to 
an energy splitting between its resonant states of about 
8meV; this allows one to "prepare" the electronic sys- 
tem in a coherent superposition of these resonant states 
such to localize the electrons within the left well |18j . 
and then to investigate the subsequent coherent charge 
and current oscillations between the two wells. More 
specifically, in order to investigate the non-trivial inter- 
play between steady-state transport properties and ul- 
trafast optical excitations, we have chosen as initial con- 
dition of the simulated experiments a single-particle den- 
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FIG. 2. (Color online) Time evolution of the photoexcited 
carrier density across a realistic GaAs/AlGaAs triple-barrier 
nanodevice (well width a w — 8 nm, barrier width a/, — 2 nm 
and height Vo = 300 meV). At t = photoexcited carriers are 
fully localized in the left well (lower panel). The subsequent 
transient dynamics (upper panel), corresponding to the solid- 
curve result in Fig. [3b, clearly shows an interwell tunneling 
dynamics characterized by a strong interplay between phase 
coherence and energy relaxation/decoherence processes (see 
text). 



3 
2 
1 

-1 





-1 

0.0 




0.5 1.0 1.5 2.0 

time (ps) 



2.5 



FIG. 3. (Color online) Ultrafast electro-optical response of a 
realistic GaAs/AlGaAs triple-barrier nanodevice (see lower 
panel in Fig. [2J) sandwiched between its electric contacts: 
current density j across the central barrier (in units of its 
Landauer-Biittiker value jo) as a function of time correspond- 
ing to an initial coherent superposition of the two resonant 
electronic states of the triple-barrier structure (ei — 48meV, 
£2 — 56 meV) and to a quasiequilibrium carrier injection from 
the left contact only (total charge density n ~ 7 x 10 16 cm~ 3 ) 
in the mesoscopic limit (a), and in the presence of inelastic 
scattering processes (b), for three different values of the effec- 
tive device length (I = 50 nm (dashed curves), I — lfim (solid 
curves), and I — 20/im (dotted curves)). 



sity matrix p aia2 given by the steady-state solution p^ s lQ , 2 
plus a contribution Ap aia2 corresponding to the coher- 
ent superposition of the resonant states just mentioned, 
which mimics the effect of an ultrafast interband excita- 
tion |18l 119] . Starting from this initial condition, we have 



investigated the value of the current density across the 
central barrier as a function of time for different values 
of the device effective lengths in the absence (mesoscopic 
limit) and presence of inelastic scattering processes, with 
varying the coupling with reservoirs. To this end, we 
have expressed such coupling in terms of an effective 
device length I, such that T r ° s — \v a \/l, where v a is 
the carrier group velocity. As shown in Fig. [3j in the 
mesoscopic limit (a) the initial current oscillations - 
unambiguous fingerprint of a coherent interwell tunnel- 
ing dynamics — are progressively suppressed due to the 
presence of the external reservoirs, and, as expected, the 
dissipation/decoherence time-scale decreases with the ef- 
fective device length I. In the presence of nonelastic scat- 
tering processes (b) the mesoscopic-limit scenario in (a) is 
partially modified: as expected, in the short-device limit 
(dashed curve) non-elastic scattering plays a minor role 
(compared to the action of the reservoirs), while in the 
opposite limit (dotted curve) the latter dominates, giving 
rise to a damping of the coherent oscillations as well as to 
a suppression of the steady-state current; in the interme- 
diate regime (solid curves) one deals with a non-trivial in- 
terplay between reservoir- and environment-induced dis- 
sipation versus decoherence, resulting in a faster damping 
of the oscillations. 

In summary, we have proposed a conceptually unified 
and physically reliable framework for the microscopic de- 
scription of energy dissipation and decoherence in open 
quantum systems, able to overcome a few crucial draw- 
backs and limitations of current quantum-device model- 
ing strategies. More specifically, we have derived a gen- 
eralized density-matrix equation, particularly suited for 
the design and optimization of new-generation semicon- 
ductor nanodevices with open spatial boundaries. The 
power and flexibility of this new modeling paradigm has 
been confirmed via fully three-dimensional simulated ex- 
periments. 
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